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A very brief status of next-to-leading order QCD calculations is given. As an example the 
next-to-leading order QCD calculations to the pp — >■ ttjj processes at the CERN Large 
Hardon Collider are presented. Results for integrated and differential cross sections are 
shown. They have been obtained in the framework of the Helac-Nlo system. 

1 Introduction 

At the CERN Large Hardon Collider we hope to uncover the mechanism of electroweak symme- 
try breaking and to find signals of physics beyond the Standard Model. Signal events have to 
be dug out from a bulk of background events, which are due to the Standard Model processes, 
mostly QCD processes accompanied by additional electroweak bosons. These constitute of final 
states with a high number of jets or identified particles. LHC data can only be meaningfully an- 
alyzed if a plethora of Standard Model background processes are theoretically under control. If 
one is content with a leading order description of the multijets final states it is possible to go to 
quite high orders, say 8-10 particles in the final states, which have to be well separated to avoid 
phase space regions where divergences become troublesome. Quite a number of tools enable one 
to do this: ALPGE^Q [1], Amegic++/Sherp41 WM, Comix/Sherp^I S], CompHEfQ [5], 
Helac-PhegasEI IMH], MadGraph/MadEventR pITU] and O'mega/WhizareIiI [IT]. Those 
tools, which are based on Feynman diagrams suffer from inefficiency at large n, where n is a 
number of particles, because the number of diagrams to be evaluated grows rapidly as n in- 
creases. Some of the tools above use methods designed to be particularly efficient at high 
multiplicities. Namely, they build up amplitudes for complex processes using off-shell recursive 
methods. Nevertheless, all these tools are completely self contained and provide amplitudes 
and integrators on their own. Although multijet observables can rather easily be modeled at 
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leading order, this description suffers several drawbacks. Leading order calculations depend 
strongly on the renormalisation scale and can therefore give only an order of magnitude esti- 
mate on absolute rates. Besides normalization, sometimes also shapes of distributions arc first 
known at higher orders. Secondly, for many scale processes like e.g. ttH, tt + nj, nj, V + nj, 
VV + nj, where V stands for and/or Z, a proper scale choice is problematic. For some 
observables dynamical scales seem to work better, for others fixed scales are applied. How 
do we know which scale to choose ? Moreover, at leading order a jet is modeled by a single 
parton, which is a very crude approximation. The situation can significantly be improved by 
including higher order corrections in perturbation theory. Next-to-leading order programs can 
be divided into three categories. Libraries with a specific list of processes at hadron-hadron 
colliders like eA). MCF]v|f| for processes with heavy quarks and/or heavy electroweak bosons. 
Nlo Jets-) — h □ [12] for jet production and Vbfnlo lij [13] for vector-boson fusion processes. 
Automatic tools based on Passarino-Vcltman [M] reduction of one-loop amplitudes for general 
2 n processes like e.g. FeynArts/FormCalc/LoopTool^ [111112] and GOLEivEl [TT] . 
And finally, automatic tools based on OPP reduction and other unitarity-based methods, which 
aim at multiparticle processes at hadron coUiders like BlackHat/Sherpa, Rocket/Mcfm 
and Helac-Nlo. Thanks to all these methods and developments several 2 — )■ 4 processes 
have recently been calculated at next-to-leading order QCD, including pp — >■ tthh [T5H20) . 
pp^W^+ 3j [HHH], pp{qq) bbbb ^,pp^ ttjj ^,pp->-Zj*+ 3j ^,pp->- W^W^jj 
production via weak-boson fusion j^l], QCD- mediated pp — W^W~^jj process [3D], and finally 
pp W+W'bb [3n[32] . Additionally, the first 2^5 process pp + 4j [33] has recently 

been calculated in the leading-color approximation. 

In this contribution, a brief report on the Helac-Nlo approach and the pp — > ttjj compu- 
tation is given. 



2 Details of the next-to-leading order calculation 

The next-to-lcading order results are obtained in the framework of Helac-Nlo based on the 
Helac-Phegas leading-order event generator for all parton level, which has, on its own, already 
been extensively used and tested in phenomenological studies see e.g. [34U38j . The integration 
over the fractions xi and X2 of the initial partons is optimized with the help of Parno |39j . 
The phase space integration is executed with the help of KaleiF^ [40] and cross checked with 
Phegas [7], both general purpose multi-channel phase space generators. The next-to-leading 
order system consists of: 

1. CutTool^^ [41], for the OPP reduction of tensor integrals with a given numerator to a 
basis of scalar functions and for the rational parts [421 - 144] : 

2. Helac-ILoop [45] for the evaluation of one loop amplitude, more specifically for the 
evaluation of the numerator functions for given loop momentum (fixed by CutTools); 
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Figure 1: Scale dependence of the total cross section for pp — > ttjj -\- X at the LHC with 
fJ'R = f-F = C ' ™t • ^^ft po-nel: The blue dotted curve corresponds to the leading order whereas 
the red solid to the next-to-leading order one. Right panel: The blue dotted curve corresponds 
to the leading order whereas the green dashed to the next-to-leading result with a jet veto of 50 
GeV. 



3. QneLQofP^ [45l|46], a library of scalar functions, which provides the actual numerical 
values of the integrals. 

4. Helac-Dipole^3 IH], automatic implementation of Catani-Seymour dipole subtraction 
[IH]: for the calculation of the real emission part. 

Let us emphasize that all parts are calculated fully numerically in a completely automatic 
manner. 



3 Numerical Results 

We consider proton-proton collisions at the LHC with a center of mass energy of y/s = 14 TeV. 
The mass of the top quark is set to be mt = 172.6 GeV. We leave it on-shell with unrestricted 
kinematics. The jets are defined by at most two partons using the kx algorithm [IHISD], with 
a separation AR = 0.8, where 

Ai?y = = l/21n[(i?, (1) 

being the rapidity and the azimuthal angle of parton i. Moreover, the recombination is only 
performed if both partons satisfy \yi\ < 5 (approximate detector bounds). We further assume 
that the jets are separated by AR = 1.0 and have lyjotl < 4.5. Their transverse momentum is 

http : //helac-phegas . web . cern . ch/helac-phegas/OneLQop . html 
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Figure 2: Distribution in the transverse momentum px of the I''* hardest jet (left panel) and the 
2"^* hardest jet (right panel) for pp — >■ ttjj + X at the LHC. The blue dotted curve corresponds 
to the leading order whereas the red solid to the next-to-leading order one. 



required to be larger than 50 GcV. We consistently use the CTEQ6 set of parton distribution 
functions [5T1[5^ . i.e. we take CTEQ6L1 PDFs with a 1-loop running ag in leading order and 
CTEQ6M PDFs with a 2-loop running at next-to-leading order. 

We begin our presentation of the final results of our analysis with a discussion of the total 
cross section. For the central value of the scale, pL^ = fip = rrit, we have obtained: 

CTLO = (120.17 ±0.08) pb, CTNLO = (106.95 ±0.17) pb. (2) 

From the above result one can read a K factor K = 0.89 which corresponds to the negative 
corrections of the order of 11%. 

The scale dependence of the corrections is illustrated in Figure [TJ We observe a dramatic 
reduction of the scale uncertainty while going from leading order to next-to-leading order. 
Varying the scale up and down by a factor 2 changes the cross section by +72% and —39% in 
the leading order case, while in the next-to-leading order case we have obtained a variation of 
— 13% and —12%. The third jet, which stems from real radiation, has not been restricted in 
this case. 

Therefore, in the next step, we study the impact of a jet veto on the third jet, which is 
simply an upper bound on the allowed transverse momentum, px- The total cross section with 
a jet veto of 50 GeV is 

CTNLo(PT,x < 50 GeV) (76.58 ± 0.17) pb , (3) 

which corresponds to K = 0.64 and negative corrections of the order of 36%. In this case a 
scale variation of —54% and —0.3% has been reached, see Figure[T]for graphical representation. 

While the size of the corrections to the total cross section is certainly interesting, it is crucial 
to study the corrections to the distributions. In Figure[2]the transverse momentum distributions 
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of the hardest and second hardest jet are shown for the pp — >■ tijj + X process. The bhie dashed 
curve corresponds to the leading order, whereas the red sohd one to the next-to-lcading order 
result. The histograms can also be turned into dynamical K-factors, which we display in the 
lower panels. Distributions demonstrate tiny corrections up to at least 200 GeV, which means 
that the size of the corrections to the cross section is transmitted to the distributions. On the 
other hand, strongly altered shapes are visible at high px especially in case of the first hardest 
jet. Let us underline here, that corrections to the high px region can only be correctly described 
by higher order calculations and arc not altered by soft-coUinear emissions simulated by parton 
showers. 

4 Summary 

We report the results of a next-to-leading order simulation of top quark pair production in 
association with two jets. With our inclusive cuts, we show that the corrections with respect 
to leading order are negative and small, reaching 11%. The error obtained by scale variation is 
of the same order. 
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